Doping and temperature dependence of electronic Raman response in cuprate 

superconductors 
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The doping and temperature dependence of the electronic Raman response in cuprate supercon- 
ductors is studied within the kinetic energy driven superconducting mechanism. It is shown that the 
temperature dependent depletion at low-energy shifts is faster in the Bi g symmetry than in the B^g 
symmetry. In analogy to the domelike shape of the doping dependent superconducting transition 
temperature, the maximal peak energy in the Bi g channel occurs around the optimal doping, and 
then decreases in both underdoped and overdoped regimes. Moreover, the overall density of Cooper 
pairs increases with increasing doping in the underdoped regime. 



Cuprate superconductors are doped Mott insulators 
with the strong electron correlation dominating the entire 
phase diagram [lj. After over 20 years extensive stud- 
ies, it has become clear that superconductivity in doped 
cupratcs results when electrons pair up into Cooper pairs 
as in the conventional superconductors How- 
ever, the superconducting (SC) transition temperature 
T c is strongly dependent on doping, and takes a dome- 
like shape with the underdoped and overdoped regimes 
on each side of the optimal doping where T c researches 
its maximum. Moreover, as a natural consequence of the 
unconventional SC mechanism 4 1 , the Cooper pairs have 
a dominant d-wave symmetry 1 , 0] . Since many physi- 
cal properties of the two-particle electron dynamics have 
been attributed to particular characteristics of low energy 
excitations determined by the electronic Raman response 
(ERR) @, i], in this Case, cl central issue to clarify the 
nature of the electron dynamics is how ERR evolves with 
doping and temperature. 

Experimentally, by virtue of systematic studies using 
the ERR measurement technique, some essential features 
of the evolution of the two-particle excitations in cuprate 
superconductors with dopingand temperature in the SC 
state have been established |5l-[la|: (a) the B\ g orienta- 
tion projects out excitations around the antinodal region 
in the Brillouin zone where the energy gap is maximal; 
(b) the Big orientation projects out excitations around 
the nodal region where the energy gap vanishes; (c) as a 
consequence of both cases (a) and (b) , ERR has a peak 
at the antinodes in the B\ g symmetry and at slightly 
lower energy at the nodes in the Bi g symmetry; how- 
ever, (d) in the underdoped regime, the peak in the B\ g 
symmetry shifts to a higher energy with decreasing dop- 
ing and tracks of the antinodal gap, while the peak in 
the Big symmetry shows the same doping dependence as 
T c and shifts to a lower energy with decreasing doping; 
and (e) the Bi g response depends linearly on energy w 
in the limit to ~ for the energy gap vanishing along 
with the diagonal directions of the Brillouin zone, while 
for the B\g orientation, the Raman vertex vanishes along 
with the energy gap at the same directions, which yields 
an additional uj 2 contribution from the line nodes of the 
vertex, and the resulting ERR varies as lu 3 . Theoreti- 
cally, an agreement has emerged that a simple d-wave 



Bardeen-Cooper-Schrieffer (BCS) formalism is useful in 
the phenomenological description of ERR in cuprate su- 
perconductors |q [l7l - r2ll |. In particular, incorporating the 
d-wave BCS formalism with the spin-fluctuation model 
[l9l - l2ll | , the role of the collective modes in the ERR ex- 
periment on cuprate superconductors is studied |2(j , and 
the result shows that the collective modes of the ampli- 
tude fluctuation of the d-wave gap yield a broad peak 
above the threshold in the B\ g spectrum. Moreover, it 
has been shown that ERR in the B\ g symmetry allows 
one to distinguish between phonon-mediated and mag- 
netically mediated d-wave superconductivity [2lj |. How- 
ever, to the best of our knowledge, ERR in cuprate su- 
perconductors has not been treated starting from a mi- 
croscopic SC theory, and no explicit calculation of the 
doping dependence of ERR has been made so far. Very 
recently, using a simple relationship between the inte- 
grated ERR and density of Cooper pairs (DOCP), it is 
shown experimentally [22J that in a similar way to the 
superfluid density [23], the overall DOCP strongly in- 
creases with increasing doping in the underdoped regime 
[22| . Thus this recently discovered doping dependence of 
the overall DOCP is also calling for an explanation. 

In this Letter, we study ERR in cuprate superconduc- 
tors based on the kinetic energy (KE) driven SC mech- 
anism [24|. We evaluate explicitly the ERR function in 
terms of the Raman density-density correlation function 
(RDDCF), and qualitatively reproduce some main fea- 
tures of ERR in cuprate superconductors 043 H3| ■ I n 
particular, we show that the temperature dependent de- 
pletion at low-energy shifts is faster in the B\ g symmetry 
than in the Bi g symmetry. Moreover, in analogy to the 
domelike shape of the doping dependence of T c , the maxi- 
mal peak energy in the Bi g channel occurs around the op- 
timal doping, and then decreases in both underdoped and 
overdoped regimes. Furthermore, the overall DOCP in- 
creases with increasing doping in the underdoped regime. 

We start from the t-J model [l|, |j|, 

H = —t C\ a C l+ f ]a + t' Cl a Ci + iy a 
if}a if]' a 

+ »j2 c l Ci « + J I2 Si - Si +v> 

ia if/ 

acting on the Hilbert subspace with no doubly occupied 
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site, i.e., G^G^ < 1, where fj = ±x, ±y, fj' = ±x±y, 
C\ a (Cia) is the creation (annihilation) operator of an 
electron with spin a, Si = (Sf,Sf,Sf) are spin opera- 
tors, and /i is the chemical potential. To deal with the 
constraint of no double occupancy, the charge-spin sepa- 
ration (CSS) fermion-spin theory (25|, [26| has been devel- 
oped, where the physics of no double occupancy is taken 
into account by representing the constrained electron op- 
erators as Cif = h\^S^ and di = h\^S^, with the spin- 
ful fermion operator h ia = e~ l ^ ia hi that describes the 
charge degree of freedom together with some effects of 
spin configuration rearrangements due to the presence of 
the doped hole itself, while the spin operator S% repre- 
sents the spin degree of freedom, then the electron local 
constraint for the single occupancy is satisfied in analyti- 
cal calculations. In this CSS fermion-spin representation, 
the t-J model Q]) is expressed as, 



H 



if} 



if}' 



ft ^ h i(T hi a + J c ff ^ Si ■ Si+jj, 



(2) 



with J ff = (1 — S) 2 J, and S = {h/ irT hi a ) — (h\hi) is the 
hole doping concentration. 

For an understanding of the physical properties of 
cuprate superconductors in the SC state, the KE driven 
SC mechanism has been developed [24j based on the 
CSS fermion-spin theory, where the interaction between 
charge carriers and spins from the KE term in the t- 
J model @ induces the d-wave charge carrier pairing 
state by exchanging spin excitations, then the electron 
Cooper pairs originating from the charge carrier pairing 
state are due to the charge-spin recombination, and their 
condensation reveals the SC ground-state. In particular, 
this SC state is a conventional BCS-like with the d-wave 
symmetry [26|,[27]], so that the basic d-wave BCS formal- 
ism is still valid in discussions of the low energy electron 
quasiparticle excitations of cuprate superconductors, al- 
though the pairing mechanism is driven by KE by ex- 
changing spin excitations. Following these previous dis- 
cussions [13, H3 , the full charge carrier Green function 
can be obtained in the Nambu representation as, 



g(k,ito n ) = Z hF 



iuj n T + £ k r 3 - A hz (k)ri 



llU)„ 



(3) 



where tq is the unit matrix, t% and T3 are Pauli ma- 
trices, the charge carrier quasiparticle spectrum Ehu = 

£|+ I A^^(k) 2 with the renormalized d-wave charge 
carrier pair gap function A^^(k) = A.f l z[cosk x — cosfe„]/2, 
while the charge carrier quasiparticle coherent weight 
ZhF and other notations are defined as same as in Ref. 
[27j , and have been determined by the self-consistent cal- 
culation [H^- 

In the CSS fermion-spin representation [25[, the elec- 
tron Green's function, 



G(k,iuj n ) 



Gii(k,Kj„), G 12 (k,iw„) 
G 2 i(k,wj„), G 22 (k, iuj n ) 



(4) 



is a convolution of the spin Green's function and charge 
carrier Green's function ([3]), and its diagonal and off- 
diagonal components Gii(i—j,t—t') = ((Ci a (t);Cj a (t'))) 
and G 2 i(i-j,t-t') = ((C^(t); Cj^(t')}) have been given 
in Ref. [13 . 

The ERR function 5(q, to) is obtained from the imag- 
inary part of RDDCF x(q, to) as [6], 



5(q,w) = [1 + ra B (w)]Imx(q,w), 



(5) 



where n B (cj) is the boson distribution function, while 
RDDCF x(q, w) is defined as, 

X(q,r-r') = -(Tp 7 (q,TK(-q,T')]), (6) 

where p 7 (q) = J^k 7kG^. + qT3Gk_a is the Raman density 
operator in the Nambu representation, with the bare Ra- 
man vertex 7k has been classified by the representations 
Big, B 2g , and A\ g of the point group as [6Hl8j. 

{b Ui ,u 3 [cos(k x a) - cos(kya)} /4, B lg , 
bLi,uj s s i n (kxa)sm(kya), B 2g , (7) 

au>i,u) s [cos(fc x a) + cos(k y a)} /4, A lg , 

respectively, where as a qualitative discussion, the mag- 
nitude of the energy dependence of the prefactors b, b' 
and a can be rescaled to units. In this case, x(q, to) can 
be obtained in terms of the electron Green's function Q 
as, 



X7172 (Qi ^) — 2 ^ 7V3 



k,p,p' 

, +l/ A fei? (p + k 



EB p B p > 
7ik72k 



q)A fez (p' + k) 

4£'/ip+k+q-E/ip' +k 



Mp> p + k + <*) V (p'> p' + k ) K\ P +k+ q ^ ';j.p'+k 

n F [(-iy'E hp , +k - (-!)»'+"' w p ,] - n F [(-l)-^ p+k+q - (-1)"+"^] 

to + (-ir £v +k - {-iy+»'u p , - (-i)^ p+k+q - {-iy+ v to v 1 1 ' 
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where the electron quasiparticle coherent weight Zp = 
jif(w) is the fermion distribution function, 
fifck = (1 + &/^«0/2 and [/ 2 \ k = (1 - &/E hk )/2, 
while the spin excitation spectrum w p and i? p have been 
given in Ref. [27|. 
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FIG. 1: (a) Bi g and (b) i?2 9 spectra as a function of energy at 
p = 0.15 for T = 0.002J. The dashed lines are a cubic and a 
linear fit for the low-energy Bi g and B^g spectra, respectively. 
Inset: the corresponding experimental results taken from Ref. 
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In cuprate superconductors, although the values of J 
and t are believed to vary somewhat from compound to 
compound 1], however, as a qualitative discussion, the 
commonly used parameters in this paper are chosen as 
t/J = 2.5 and t'/t = 0.3. In this case, we have per- 
formed a calculation for the ERR function ([3]) in both 
Big and B 2g orientations, and the results of (a) the B\ g 
and (b) B 2g spectra with doping p = 0.15 at tempera- 
ture T — 0.002J are plotted in Fig. Q] For comparison, 
the corresponding experimental results [|| are also pre- 
sented in Fig. Q] (inset). Obviously, both B\ g and B 2g 
spectra are characterized by the presence of the pair- 
breaking peaks, however, the peak in the B\ g geometry 
develops at roughly 30% higher energy than that in the 
B 2g geometry. Moreover, we have also fitted our present 
results, and found that the low-energy spectra almost 
rise as w 3 in the B\ g channel and linearly with uj in the 
B 2g channel. However, in contrast to the conventional 
superconductors, the intensities of ERR in both B\ g and 
B 2g channels are equal to zero when energy uj — 0, in 
other words, there are no sharp onset of the intensities 
at a threshold. This can be understood from the physical 



property of RDDCF ©. In low temperatures, the spins 
center around the [±7r, ±7r] points, then the main contri- 
bution from the spins comes from the [±7r, ±7t] points. 
In this case, RDDCF ([SJ in the momentum transfers 
q ~ limit can be approximately reduced in terms of 
w p=[±tt ±tt] ~ and one of the self-consistent equations 
1/2 



(S+Sr) = (l/A0£ p £ p coth(/3 Wp /2)/(2 Wp ) as 



1 A 2 (k) 1 

X 7 i7 2 (q~ 0,w) « Zp— 2^7ik72k z e2 tanh[-/3£ k ] 



1 1 

+ 



2£ k + oj 2£ k - uj J ' 



(9) 



where the renormalized d-wave electron pair gap func- 
tion Az(k) w A/jz(k + kA) with kA = [7r,7r], and the 
electron quasiparticle spectrum i? k w i?h k + kA ■ Since the 
d-wave gap function vanishes along the diagonal direc- 
tions of the Brillouin zone, then as seen from Eq. 
this leads to an absence of the threshold in cuprate su- 
perconductors. These results are qualitatively consistent 
with the experimental results [7H9|]- 
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FIG. 2: (a) B\ g and (b) Big spectra as a function of energy 
at p = 0.20 for T = 0.03 J (solid line), T = 0.05 J (dotted 
line), and T — 0.07J (dashed line). Inset: the corresponding 
experimental results taken from Ref. [l6l . 

For a better understanding of the evolution of ERR 
with temperature, we have further performed a calcula- 
tion for the ERR function ([5]) with different tempera- 
tures, and the results of 5(q, w) in (a) the B\ g and (b) 
B 2g channels as a function of energy with T = 0.03 J 
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FIG. 3: (a) Bi g and (b) Big spectra as a function of energy at 
p = 0.20 (solid line) andp = 0.22 (dotted line) for T = 0.002 J. 
Inset: the corresponding experimental results taken from Ref. 
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(solid line), T = 0.05 J (dotted line), and T = 0.07 J 
(dashed line) for p = 0.20 are plotted in Fig. [2] in com- 
parison with the corresponding experimental results [l6l | 
(inset). Our results show that there is a strong tempera- 
ture dependence of the peaks in both B\ g and Big orien- 
tations in the overdoped regime, where both B\ g and Bi g 
peaks at low temperatures soften in energy and decreases 
in intensity as the temperature is raised. Moreover, as 
seen from Eq. ©, the intensity of the peak follows a 
pair gap type temperature dependence, and disappears 
at T c . Simultaneously, the low energy continuum grows 
in agreement with the transfer of spectral weight. Fur- 
thermore, the temperature dependent depletion at low- 
energy shifts is faster in the B\ g symmetry than in the 
Bi g symmetry. Although the strong temperature depen- 
dence of the Big spectrum throughout the SC dome and 
Big spectrum in the overdoped regime are qualitatively 
consistent with the corresponding experimental results 
[T5I Il6j |. the temperature dependence of the B\ g spec- 
trum in the underdoped regime is in disagreement with 
the corresponding experimental result [l5l Il6j , where the 
nontemperature dependence (or weakly increases as the 
temperature is raised) is observed. This reflects that the 
excitations observed in the B\ g channel is only related 
to superconductivity in the overdoped regime, and then 
ERR in the B\ g and Big channels become drastically dis- 
tinct in the underdoped regime, which can be seen more 
clearly in the following discussions of the doping depen- 



dence of ERR in the B\ g channel. 
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FIG. 4: (a) B\ g and (b) Big peaks as a function of doping for 
T — 0.002J. Inset: the corresponding experimental results 
taken from Ref. [la. 
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FIG. 5: Doping dependence of (a) the density of Cooper pairs 
in the B\ g (solid line) and Big (dashed line) channels, and (b) 
the sum of the densities of Cooper pairs in the B\ g and Big 
channels for T = 0.002J. Inset: the corresponding experi- 
mental results taken from Ref. 1 2 21 . 
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To analyze the evolution of ERR in Fig. [T] with dop- 
ing, we have made a series of calculations for the ERR 
function ([5]) with different dopings, and the results of 
(a) the Big and (b) Bi g spectra as a function of en- 
ergy with p = 0.20 (solid line) and p = 0.22 (dotted 
line) at T = 0.002 J are plotted in Fig. [3] in comparison 
with the corresponding experimental results [l2[ (inset). 
It is shown clearly that although the intensities of the 
peaks in both B\ g and B<iq symmetries are strongly sup- 
pressed in the overdoped regime when reducing doping, 
the peak energies change drastically, i.e., the peak ener- 
gies in both B\g and Bi g channels decrease with increas- 
ing doping. Furthermore, we have found that the peak 
energies in both B\ g and Big channels increase with in- 
creasing doping in the underdoped regime. To show this 
point clearly, we have calculated the ERR function ([5]) 
throughout the SC dome, and then employed the shift 
of the leading-edge mid-point as a measurement of the 
magnitude of the pair gap at each doping just as it has 
been done in the ERR experiments [lOl - Rq ]. The results 
for the extracted (a) B\ g and (b) Bi g peaks as a func- 
tion of doping with T = 0.002J are plotted in Fig. g] 
For comparison, the corresponding experimental results 
ilSj are also presented in Fig. [4] (inset). The result of 
the peak energy in the Bi g symmetry continuously fol- 
lows T c throughout the SC dome as oj^^ oc T c , i.e., the 

maximal peak energy w^ o 2 ^ k in the Big symmetry occurs 
around the optimal doping, and then decreases in both 
underdoped and overdoped regimes. In this sense, it can 
be identified with the SC gap in the whole doping range, 
in qualitative agreement with the experimental results 
[10l - [la |. However, in the B\ g orientation, although the 
doping dependence of the peak energy is qualitatively 
consistent with the experimental data in the overdoped 
regime, it is in disagreement with the experimental data 
in the underdoped regime [icil - fl5l |. where the peak en- 
ergy increases essentially linearly as doping is reduced. 
This reflects that the peak energy in the B\ g symmetry 
is progressively disconnected from superconductivity as 
one goes from the overdoped regime to the underdoped 
regime. 

Now we turn to discuss the evolution of DOCP with 
doping. In the case of the low temperatures and momen- 
tum transfers q ~ limit, the integral of the imaginary 
part of RDDCF over energy oj can be obtained from Eq. 
© as, 

Ppah-s = / dwx 7172 (q~0,u;) 
Jo 

« 27rZ|l^ 7lk72k ^^tanh[i/3£ k ]. (10) 

However, the sum $Z k (A|(k) / '£7 k )tanh(/3.E k /2) is equal 
to 4X; k (C4Vk) 2 tanh(/3£; k /2), where the electron quasi- 
particle coherence factors U£ ~ ^2/ik+k A an< ^ ^k ~ 
Uf hk+k describe the probabilities of the Cooper pair be- 
ing occupied and unoccupied, respectively. Furthermore, 



this sum is not vanishing only around the Fermi energy 
in the rang e of 2Az(k). This quantity corresponds to 
DOCP [22], formed around the Fermi surface as the en- 
ergy gap is opening. In this sense, /3 pa irs is proportional to 
DOCP, weighted by the Raman vertex [22|]. In this case, 
we have performed a calculation for the doping depen- 
dence of Ppairs, and the results of (a) Ppairs ( sou d line) 
and Ppafrs (dashed line) and (b) the sum of pp*f IS and 

Ppafrs as a function of doping with T = 0.002 J are plotted 
in Fig. [5] in comparison with the corresponding experi- 
mental data [22j (inset). The present results show that 
in the momentum space, p pa i rs is strongly anisotropic 

as a function of doping. Although /Op a f rs is much larger 

than Ppaf rs in the corresponding doping, they almost lin- 
early increase with increasing doping in the underdoped 
regime. In particular, at low doping, p^ rs and Ppff rs be- 
come very small and vanishc below p ~ 0.05. Moreover, 
in a similar way to the superfluid density, the sum of 
P^fr= a nd p^? , which describes the overall DOCP, in- 
creases linearly with increasing doping in the underdoped 
regime, in qualitative agreement with the experimental 
data [221 . 

A nature question is why the doping and temperature 
dependence of the B± g spectrum in the overdoped regime 
and Big spectrum throughout the SC dome in cuprate 
superconductors can be described qualitatively within 
the KE driven SC mechanism. The reason is that the 
KE driven SC state is the conventional d-wave BCS like 
[13, H3- I n particular, this KE driven d-wave SC state 
is controlled by both SC gap and quasiparticle coher- 
ence, then the maximal T c (then gap parameter) occurs 
around the optimal doping, and decreases in both under- 
doped and overdoped regimes. On the other hand, based 
on the phenomenological BCS formalism with a simple 
d-wave gap A k = A[cosk x — cosk y ], it has been shown Q 
that the positions of the maxima in the B\ g and Big scat- 
tering both scale with the gap parameter A. However, 
in the present paper, the ERR function of cuprate su- 
perconductors ([5]) is obtained in terms of the imaginary 
part of RDDCF © within the KE driven SC mecha- 
nism, although its expression is similar to that obtained 
within the phenomenological d-wave BCS formalism Q. 
This is why the doping and temperature dependence of 
the Big spectrum in the overdoped regime and Big spec- 
trum throughout the SC dome in cuprate superconduc- 
tors can be described qualitatively within the KE driven 
SC mechanism. 

Although many different theories have been put for- 
ward in order to explain the nature of ERR in the B\ g 
orientation in the underdoped regime, however, its full 
understanding is still controversial. One point of view is 
that the experimental result from ERR [HHLa ] does not 
necessary imply two distinct gaps, but still it does imply 
that the gap behaves differently in the nodal and antin- 
odal momentum regimes. In particular, it has been ar- 
gued that the ERR data in both B\ g and Bi g geometries 
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in the underdoped regime can be equally well explained 
within a one-gap scenario if final-state interactions are 
taken into account [28]. In other theories, the normal 
state pseudogap is associated with the nature of the B\ g 
peak energy in the underdoped regime [29j. However, in 
the present case, the continuation of the treatments from 
the SC to normal states is not straightforward. This can 
be seen directly from Eq. ©, the intensity of ERR van- 
ishes proportional to A|(k) as T T c . In this case, 
ERR in the underdoped regime has been studied within 
the fluctuation-exchange approximation for the Hubbard 
model by phenomenologically introducing a pseudogap 
(30| , where this pseudogap leads to that a corresponding 
peak evolves continuously in the B\ g spectrum as T de- 
creases in the normal state and below T c . However, this 
theory can not be applied to discuss the doping depen- 
dent behaviors of ERR in B\ g and £?2 5 symmetries. 

Within the KE driven SC mechanism, our present re- 
sult of the doping and temperature dependence of ERR 
in the B\ g symmetry is in disagreement with the exper- 
imental data 0-QJ, Hl| in the underdoped regime. In 
Ref. [3, it has been argued that for a correct descrip- 
tion of the Bi g spectrum in the underdoped regime, two 
essential ingredients should be taken into account within 
the d-wave BCS formalism: (a) the quasiparticle spec- 
tral weight Zf(\s.) as well as the vertex correction; and 
(b) a general gap also should be taken into account by 
including the higher order of the harmonic component in 
the simple d-wave gap function Ak = A[cosfc K — cosk y ]. 
In our present studies, the first condition is partially sat- 
isfied, since the doping dependence of the quasiparticle 
spectral weight Zp has been included in our above dis- 
cussions. However, the vertex correction for the inter- 
action of the quasiparticles due to the presence of the 
spin fluctuation is not included, which also is doping and 
momentum dependent. For the second condition, a more 
general gap Ak = A[cos/c x — cosfcj, + B(cos2k x — cos2fc a )] 



has been observed [3l|. Although this gap function still 
is basically consistent with the d-wave symmetry, it is 
obviously that there is a significant deviation from the 
simple d-wave form. On the other hand, the final-state 
interactions are also very important for interpreting the 
unusual Bi g _spectrum in the underdoped regime as point 
out in Ref. l28l and should be taken into account. Thus 
an important issue is how these essential ingredients that 
have been dropped in the present studies are taken into 
account within the KE driven SC mechanism for a cor- 
rect description of the B\ g spectrum in the underdoped 
regime. These and the related issues are under investi- 
gation now. 

In summary, we have studied the doping and temper- 
ature dependence of ERR in cuprate superconductors 
based on the KE driven SC mechanism. Our results show 
that the temperature dependent depletion at low-energy 
shifts is faster in the B\ g symmetry than in the Bi g sym- 
metry. When increasing doping, the intensities of the 
peaks in both B\ g and Bi g orientations are strongly in- 
creased. In particular, in analogy to the domelike shape 
of the doping dependence of T c , the maximal peak energy 
in the Bi g channel occurs around the optimal doping, 
and then decreases in both underdoped and overdoped 
regimes. Moreover, in a similar way to the supcrfluid 
density, the overall DOCP increases with increasing dop- 
ing in the underdoped regime. 
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